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Abstract 

Positive results of dark matter searches in experiments DAMA/Nal 
and DAMA/LIBRA taken together with negative results of other 
groups can imply nontrivial particle physics solutions for cosmologi- 
cal dark matter. Stable particles with charge -2 bind with primordial 
helium in 0-hclium "atoms" (OHe), representing a specific Warmer 
than Cold nuclear-interacting form of dark matter. Slowed down in 
the terrestrial matter, OHe is elusive for direct methods of under- 
ground Dark matter detection like those used in CDMS experiment, 
but its low energy binding with nuclei can lead to annual variations 
of energy release in the interval of energy 2-6 keV in DAMA/Nal and 
DAMA/LIBRA experiments. Schrodinger equation for system of nu- 
cleus and OHe is considered and reduced to an equation of relative 
motion in a spherically symmetrical potential, formed by the Yukawa 
tail of nuclear scalar isoscalar attraction potential, acting on He be- 
yond the nucleus, and dipole Coulomb repulsion between the nucleus 
and OHe at distances from the nuclear surface, smaller than the size 
of OHe. The values of coupling strength and mass of meson, medi- 
ating scalar isoscalar nuclear potential, are rather uncertain. Within 
these uncertainties and in the approximation of rectangular poten- 
tial wells we find a range of these parameters, at which the sodium 
and/or iodine nuclei have a few keV binding energy with OHe. At 
nuclear parameters, reproducing DAMA results, the energy release 
predicted for detectors with chemical content other than Nal differ in 
the most cases from the one in DAMA detector. In particular, it is 
shown that in the case of CDMS germanium state has binding energy 
with OHe beyond the range of 2-6 keV and its formation should not 
lead to ionization in the energy range of DAMA signal. Due to dipole 
Coulomb barrier, transitions to more energetic levels of Na(I)-|-OHe 
system with much higher energy release are suppressed in the cor- 
respondence with the results of DAMA experiments. The proposed 
explanation inevitably leads to prediction of abundance of anomalous 
Na and I, corresponding to the signal, observed by DAMA. 
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1 Introduction 



The widely shared behef is that the dark matter, corresponding to 25% of 
the total cosmological density, is nonbaryonic and consists of new stable 
particles. One can formulate the set of conditions under which new parti- 
cles can be considered as candidates to dark matter (see e.g. [H [3] for 
review and reference): they should be stable, saturate the measured dark 
matter density and decouple from plasma and radiation at least before the 
beginning of matter dominated stage. The easiest way to satisfy these con- 
ditions is to involve neutral weakly interacting particles. However it is not 
the only particle physics solution for the dark matter problem. In the com- 
posite dark matter scenarios new stable particles can have electric charge, 
but escape experimental discovery, because they are hidden in atom-like 
states maintaining dark matter of the modern Universe. 

It offers new solutions for the physical nature of the cosmological dark 
matter. The main problem for these solutions is to suppress the abundance 
of positively charged species bound with ordinary electrons, which behave 
as anomalous isotopes of hydrogen or helium. This problem is unresolvable, 
if the model predicts stable particles with charge -1, as it is the case for 
tera-electrons [11[5]. To avoid anomalous isotopes overproduction, stable 
particles with charge -1 should be absent, so that stable negatively charged 
particles should have charge -2 only. 

Elementary particle frames for heavy stable -2 charged species are pro- 
vided by: (a) stable " antibaryons" UUU formed by anti-C/ quark of fourth 
generation P, [71 EJ [9] (b) AC-leptons (HUTOllIl], predicted in the extension 
[To] of standard model, based on the approach of almost-commutative ge- 
ometry [12j . (c) Technileptons and anti-technibaryons [13] in the framework 
of walking technicolor models (WTC) [14J. (d) Finally, stable charged clus- 
ters of (anti)quarks U5 of 5th family can follow from the approach, 
unifying spins and charges [l5] . 

In the asymmetric case, corresponding to excess of -2 charge species, 
X , as it was assumed for (UUU) in the model of stable [/-quark of 
a 4th generation, as well as can take place for (U5U5U5) in the approach 
[15] their positively charged partners effectively annihilate in the early Uni- 
verse. Such an asymmetric case was realized in [13] in the framework of 
WTC, where it was possible to find a relationship between the excess of neg- 
atively charged anti-techni-baryons (UU) and/or technileptons ^ and 
the baryon asymmetry of the Universe. The relationship between baryon 
asymmetry and excess of -2 charge stable species is supported by sphaleron 
transitions at high temperatures and can be realized in all the models, 
in which new stable species belong to non-trivial representations of elec- 
troweak SU(2) group. 
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After it is formed in the Standard Big Bang Nucleosynthesis (SBBN), 
'^He screens the X charged particles in composite {'^He'^^X ) O- 
helium "atoms" l6|. For different models of X these "atoms" are also 
called ANO-helium [71, ^8], Ole-hehum [SKII] or techni-O-hehum [13]. We'h 
call them all 0-helium {OHe) in our further discussion, which follows the 
guidelines of [16] , 

In all these forms of 0-helium X behave either as leptons or as spe- 
cific "heavy quark clusters" with strongly suppressed hadronic interaction. 
Therefore 0-helium interaction with matter is determined by nuclear in- 
teraction of He. These neutral primordial nuclear interacting objects con- 
tribute to the modern dark matter density and play the role of a nontrivial 
form of strongly interacting dark matter [17^ [T8] . The active influence of 
this type of dark matter on nuclear transformations seems to be incompat- 
ible with the expected dark matter properties. However, it turns out that 
the considered scenario of nuclear-interacting 0-helium Warmer than Cold 
Dark Matter is not easily ruled out [6l fTTj [T3l [T9] and challenges the exper- 
imental search for various forms of 0-helium and its charged constituents. 

Here we concentrate on its effects in underground detectors. We present 
qualitative conflrmation of the earlier guess [16^ |20] that the positive re- 
sults of dark matter searches in DAMA/Nal (see for review [21]) and 
DAMA/LIBRA [22j experiments can be explained by 0-helium, resolving 
the controversy between these results and negative results of other experi- 
mental groups. 

2 OHe in the terrestrial matter 

The evident consequence of the 0-helium dark matter is its inevitable pres- 
ence in the terrestrial matter, which appears opaque to 0-helium and stores 
all its in-falling flux. 

After they fall down terrestrial surface, the in-falling OHe particles are 
effectively slowed down due to elastic collisions with matter. Then they 
drift, sinking down towards the center of the Earth with velocity 

y = ^ ^SOSg^^/^cm/s. (1) 
nav 

Here yl 30 is the average atomic weight in terrestrial surface matter, 
n = 2.4 • 1^"^^ /A cm~^ is the number density of terrestrial atomic nuclei, av 
is the rate of nuclear collisions, nio ~ Mx + ^rnp = TeV is the mass of 
0-helium, Mx is the mass of the X component of 0-helium, nip is the 
mass of proton and g = 980 cm/ s'^. 

Near the Earth's surface, the 0-helium abundance is determined by the 
equilibrium between the in-falling and down-drifting fluxes. 
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The in-falling 0-helium flux from dark matter halo is 



Sir 

where Vh-speed of Solar System (220 km/s), VE-speed of Earth (29.5 km/s) 
and no = 3-10~^S'3"^ cm~'^ is the local density of 0-helium dark matter. For 
qualitative estimation we don't take into account here velocity dispersion 
and distribution of particles in the incoming flux that can lead to significant 
effect. 

At a depth L below the Earth's surface, the drift timescale is tdr ~ L/V, 
where V ~ 4005*3 cm/ s is given by Eq. ([1]). It means that the change of the 
incoming flux, caused by the motion of the Earth along its orbit, should lead 
at the depth L ~ 10^ cm to the corresponding change in the equilibrium 
underground concentration of OHe on the timescale tdr ~ 2.5 • lO'^S^^ s. 

The equilibrium concentration, which is established in the matter of 
underground detectors at this timescale, is given by 

27r • F nav ,— , 
rioE = =no——-\Vh + VE\, (2) 

V 4g 

where, with account for > Ve, relative velocity can be expressed as 



Vo\ = ^{Vh + VeY = JV^ + + VhVEsin{e) 



Here 9 = Lo{t — to) with to = 2tt/T, T = lyr and to is the phase. Then the 
concentration takes the form 

noE = n^^i + ■ sin{Lo{t - to)) (3) 

So, there are two parts of the signal: constant and annual modulation, 
as it is expected in the strategy of dark matter search in DAMA experiment 
[22] . 

Such neutral {'^He'^^X ) "atoms" may provide a catalysis of cold 
nuclear reactions in ordinary matter (much more effectively than muon 
catalysis). This effect needs a special and thorough investigation. On 
the other hand, X capture by nuclei, heavier than helium, can lead to 
production of anomalous isotopes, but the arguments, presented in [HI [TTl 
[T3] indicate that their abundance should be below the experimental upper 
limits. 

It should be noted that the nuclear cross section of the 0-helium interac- 
tion with matter escapes the severe constraints [T8j on strongly interacting 
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dark matter particles (SIMPs) [171 H] imposed by the XQC experiment 
[23j . Therefore, a special strategy of direct 0-helium search is needed, as 
it was proposed in [23] . 

In underground detectors, OHe "atoms" are slowed down to thermal 
energies and give rise to energy transfer ~ 2.5 • lO^'^eVA/S'a, far below 
the threshold for direct dark matter detection. It makes this form of dark 
matter insensitive to the severe CDMS constraints [25j. However, OHe 
induced processes in the matter of underground detectors can result in 
observable effects. 

3 Low energy bound state of O-helium with nuclei 

In the essence, our explanation of the results of experiments DAMA/Nal 
and DAMA/LIBRA is based on the idea that OHe, slowed down in the 
terrestrial matter and present in the matter of DAM A detectors, can form 
a few keV bound state with nucleus, in which OHe is situated beyond the 
nucleus. Formation of such bound state leads to the corresponding energy 
release and ionization signal, detected in DAMA experiments. 

3.1 Low energy bound state of O-helium with nuclei 

We assume the following picture: at the distances larger, than its size, OHe 
is neutral and it feels only Yukawa exponential tail of nuclear attraction, 
due to scalar-isoscalar nuclear potential. It should be noted that scalar- 
isoscalar nature of He nucleus excludes its nuclear interaction due to vr or 
p meson exchange, so that the main role in its nuclear interaction outside 
the nucleus plays a meson exchange, on which nuclear physics data are 
not very definite. When the distance from the surface of nucleus becomes 
smaller than the size of OHe, the mutual attraction of nucleus and OHe 
is changed by dipole Coulomb repulsion. Inside the nucleus strong nuclear 
attraction takes place. In the result the spherically symmetric potential 
appears, given by 



Here Afje = 4, Zue = 2 are atomic weight and charge of helium, A and Z 
are respectively atomic weight and charge of nucleus, // and are the mass 
and coupling of scalar-isoscalar meson - mediator of nuclear attraction, Tq is 
the size of OHe and F{r) is its electromagnetic formfactor, which strongly 
suppresses the strength of dipole electromagnetic interaction outside the 



U = 



AHeAg'^exp{-nr) ZmZf? 



• F{r) 



(4) 



r 



OHe 



atom" . 



5 



Schrodinger equation for this system is reduced (taking apart the equa- 
tion for the center of mass) to the equation of relative motion for the reduced 
mass 

Anipnio 

m = — , (5) 

Anip + nio 

where rrip is the mass of proton. 

In the case of orbital momentum /=0 the wave functions depend only 
on r. 

To simplify the solution of Schrodinger equation we approximate the 
potential ([4j) by a rectangular potential that consists of a deep potential well 
within the radius of nucleus Ra, of a rectangular dipole Coulomb potential 
barrier outside its surface up to the radial layer a = Ra + ro, where it is 
suppressed by the OHe atom formfactor, and of the outer potential well of 
the width ~ l//i, formed by the tail of Yukawa nuclear interaction. It leads 
to the approximate potential, given by 

r<RA:U = U^ = -^'"''P^-^''^^ 



Ra 

j-RA+To 2Z( 

Ra < r < a : U = U-2 



rRA+To 2Za4TT(ro/x) 

JRa X "-^ 



(6) 



Ra + ro 
b<r:U = U4 = 0, 

presented on Fig. [TJ 

Solutions of Schrodinger equation for each of the four regions, indicated 
on Fig. [H are considered in Appendix. In the result of their sewing one 
obtains the condition for the existence of a low-energy level in OHe-nucleus 
system, 

where A;3 and 6 are, respectively, the wave number and phase of the wave 
function in the region III (see Appendix for details). 

With the use of the potential ([6]) in the Eq.(I7|), intersection of the two 
lines gives graphical solution presented on Fig. [2j 

Based on this solution one obtains from Eq. p2p the energy levels of a 
bound state in the considered potential well. 

The energy of this bound state and its existence strongly depend on 
the parameters fi and of nuclear potential (jH). On the Fig. [3] the region 
of these parameters, giving 2-6 keV energy level in OHe bound states with 
sodium and iodine are presented. In these calculations the mass of OHe 
was taken equal to nio = ITeV. 



6 




7 




Figure 3: The region of parameters fi and , for which Na and I have a 
level in the interval 2-6 keV. For each nucleus two narrow strips determine 
the region of parameters, at which the bound system of this element with 
OHe has a level in 2-6 keV energy range. The outer line of strip corresponds 
to the level of 6 keV and the internal line to the level of 2 keV. The region 
of intersection of strips correspond to existence of 2-6 keV levels in both 
OHe-Na and OHe-I systems, while the piece of strip between strips of other 
nucleus corresponds to the case, when OHe bound state with this nucleus 
has 2-6 keV level, while the binding energy of OHe with the other nuclei is 
less than 2 keV by absolute value. 
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The rate of radiative capture of OHe by nuclei should be accurately 
calculated with the use of exact form of wave functions, obtained for the 
OHe-nucleus bound state. This work is now in progress. One can use 
the analogy with the radiative capture of neutron by proton, considered in 
textbooks (see e.g. [26]) with the following corrections: 

• There is only El transition in the case of OHe capture. 

• The reduced masses of n-p and OHe-nucleus systems are different 

• The existence of dipole Coulomb barrier leads to a suppression of the 
cross section of OHe radiative capture. 

With the account for these effects our first estimations give the rate of OHe 
radiative capture, reproducing the level of signal, detected by DAMA. 

Formation of OHe-nucleus bound system leads to energy release of its 
binding energy, detected as ionization signal in DAMA experiment. In the 
context of our approach the existence of annual modulations of this signal 
in the range 2-6 keV and absence of such effect at energies above 6 keV 
means that binding energy of Na-OHe and I-OHe systems should not exceed 
6 keV, being in the range 2-6 keV for at least one of these elements. These 
conditions were taken into account for determination of nuclear parameters, 
at which the result of DAMA can be reproduced. At these values of /u and 

energy of OHe binding with other nuclei can strongly differ from 2-6 
keV. In particular, energy release at the formation of OHe bound state 
with thallium can be larger than 6 keV. However, taking into account that 
thallium content in DAMA detector is 3 orders of magnitude smaller, than 
Nal, such signal is to be below the experimental errors. 

It should be noted that the results of DAMA experiment exhibit also 
absence of annual modulations at the energy of MeV-tens MeV. Energy 
release in this range should take place, if OHe-nucleus system comes to the 
deep level inside the nucleus (in the region I of Fig. [T]). This transition 
implies tunneling through dipole Coulomb barrier and is suppressed below 
the experimental limits. 

3.2 Energy levels in other nuclei 

For the chosen range of nuclear parameters, reproducing the results of 
DAMA/Nal and DAMA/LIBRA, we can calculate the binding energy of 
OHe-nucleus states in nuclei, corresponding to chemical composition of set- 
ups in other experiments. The results of such calculation for germanium, 
corresponding to the detector of CDMS experiment, are presented on Fig. 
m For all the parameters, reproducing results of DAMA experiment the 
predicted energy level of OHe-germanium bound state is beyond the range 
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Figure 4: Energy levels in OHe bound system with germanium. The range 
of energies close to energy release in DAMA experiment is blown up to 
demonstrate that even in this range there is no formal intersection with 
DAMA results. 



2-6 keV, being dominantly in the range of tens - few-tens keV by absolute 
value. It makes elusive a possibility to test DAMA results by search for 
ionization signal in the same range 2-6 keV in other set-ups with content 
that differs from Na and I. In particular, our approach naturally predicts 
absence of ionization signal in the range 2-6 keV in accordance with the 
recent results of CDMS [27]. 

We have also calculated the energies of bound states of OHe with xenon 
(Fig. [5]), argon (Fig. [6|), carbon (Fig. [7]), aluminium (Fig. [8]), fluorine (Fig. 
|9]), chlorine (Fig. [10]) and oxygen (Fig. [TTI) . 

3.3 Superheavy OHe 

In view of possible applications for the approach, unifying spins and charges 
[15j . we consider here the case of superheavy OHe, since the candidate for 
X , coming from stable 5th generation (u^u^u^) is probably much heavier, 
than 1 TeV. With the growth of the mass of 0-helium the reduced mass 
([5|) slightly grows, approaching with higher accuracy the mass of nucleus. 
It extends a bit the range of nuclear parameters ^ and g'^, at which the 
binding energy of OHe with sodium and/or iodine is within the range 2- 
6 keV (see Fig. [T2|) . At these parameters the binding energy of 0-helium 
with germanium and xenon are presented on figures [13] and [TU respectively. 
Qualitatively, these predictions are similar to the case of Ss = 1. Though 
there appears a narrow window with OHe-Ge binding energy, below 6 keV 
for the dominant range of parameters energy release in CDMS is predicted 
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Figure 5: Energy levels in OHe bound system with xenon. 
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Figure 6: Energy levels in OHe bound system with argon. 
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Figure 7: Energy levels in OHe bound system with carbon. 
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Figure 8: Energy levels in OHe bound system with aluminium. 
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Figure 9: Energy levels in OHe bound system with fluorine. 
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Figure 10: Energy levels in OHe bound system with chlorine. 
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Figure 11: Energy levels in OHe bound system with oxygen. 




Figure 12: The range of parameters /x and g'^, for which Na and I have a 
level in the interval 2-6 keV for 5*3 = 100. This range becomes a bit wider 
as compared with the case of S3 = 1, presented on Fig. [3l 
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Figure 13: Energy levels in OHe bound system with germanium. 



to be of the order of few tens keV. 

4 Conclusions 

To conclude, the results of dark matter search in experiments DAMA/Nal 
and DAMA/LIBRA can be explained in the framework of composite dark 
matter scenario without contradiction with negative results of other groups. 
This scenario can be realized in different frameworks, in particular in Min- 
imal Walking Technicolor model or in the approach unifying spin and 
charges and contains distinct features, by which the present explanation 
can be distinguished from other recent approaches to this problem |28j (see 
also review and more references in |29]). 

Our explanation is based on the mechanism of low energy binding of 
OHe with nuclei. We have found that within the uncertainty of nuclear 
physics parameters there exists a range at which OHe binding energy with 
sodium and/or iodine is in the interval 2-6 keV. Radiative capture of OHe 
to this bound state leads to the corresponding energy release observed as 
an ionization signal in DAMA detector. 

OHe concentration in the matter of underground detectors is deter- 
mined by the equilibrium between the incoming cosmic flux of OHe and 
diffusion towards the center of Earth. It is rapidly adjusted and follows the 
change in this flux with the relaxation time of few minutes. Therefore the 
rate of radiative capture of OHe should experience annual modulations re- 
flected in annual modulations of the ionization signal from these reactions. 

An inevitable consequence of the proposed explanation is appearance 
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Figure 14: Energy levels in OHe bound system with xenon. 



in the matter of DAMA/Nal or DAM A/LIBRA detector anomalous su- 
perheavy isotopes of sodium and/or iodine, having the mass roughly by 
nio larger, than ordinary isotopes of these elements. If the atoms of these 
anomalous isotopes are not completely ionized, their mobility is determined 
by atomic cross sections and becomes about 9 orders of magnitude smaller, 
than for 0-helium. It provides their conservation in the matter of de- 
tector. Therefore mass-spectroscopic analysis of this matter can provide 
additional test for the 0-helium nature of DAMA signal. Methods of such 
analysis should take into account the fragile nature of OHe-Na (and/or 
OHe-I) bound states. Their binding energy is only few keV. 

With the account for high sensitivity of our results to the values of 
uncertain nuclear parameters and for the approximations, made in our cal- 
culations, the presented results can be considered only as an illustration of 
the possibility to explain puzzles of dark matter search in the framework 
of composite dark matter scenario. However, even at the present level of 
our studies we can make a conclusion that the ionization signal expected 
in detectors with the content, different from Nal, can be dominantly in the 
energy range beyond 2-6 keV. Therefore test of results of DAMA/Nal and 
DAMA/LIBRA experiments by other experimental groups can become a 
very nontrivial task. 
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Appendix. Solution of Schrodinger equation for 
rectangular well 

In the 4 regions, indicated on Fig. [H Schrodinger equation has the form 

1 ^2 

^ ■ rd^^'''^^^ + hirf^'i = 0,ki{r) =h = y/2m{Ui - \E\); (8) 



1 £ 

7d^(^^2) + A;2(r)V2 = 0, k2{r) = k2 = ^/2m{U2 - \E\y, (9) 



1 d 



•2 



III : -^(^^3) + k3irfij3 = 0, ksir) = k^ = ^2m{U3 - \E\)- (10) 
1 <f 



IV : —r^{riJi) - k^irf^^ = 0, ki{r) = k^ = V2m|S|. (11) 

The wave functions in these regions with the account for the boundary 
conditions have the form [30] 

I:i., = A'^^; (12) 
r 

Bi ■ exp{-k2r) + B2 ■ exp{k2r) 
II ■■ V'2 = ; (13) 

r 

III : V's = C^^^^^^^ (14) 
r 

ly : ^4 = Z^^^PizM (15) 
r 

The conditions of continuity of a logarithmic derivative ^ = jt/^ 
at the boundaries of these regions r = i?^, r = a and r = 6 are given by 



T TT u wi. D ^ I, exp{k2RA) - F ■ exp{-k2RA) 
I-II:k,. ctg{k,RA) = k2 ■ (16) 



rr rrr i ^ (u ^ ^\ i, 6X^(^20) - F • exp{-k2a) 

II - III : ks ■ ctg{k3a + 6) = k2 --. — 1—. — r, (17) 

exp[k2a) + I ■ exp[—K2a) 

III -IV :k3- ctgik^b + 6) = -k^, (18) 
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where 

F = Bi/B2. (19) 

Now we can solve this system of equations for 3 variables. It follows 
from Eq. (jl6p that 

T? ^o; D ^ ^2 - A;i • ctg{kiRA) , . 

F = exp{2k2RA) ■ T—T , p V 

k2 + ki ■ ctg[kiRA) 

and from Eq. (fT7|) 

S = -ksa + arcctg{— ■ ^^^(^20) - F ■ exp{-k2a) . .^x) 
ks exp{k2a) + F ■ exp{—k2a) 



Since 



E = mr)-^, (22) 

one has 

ki = ^2mU -kl, (23) 
Then Eq. ijlSp has the form 

[ ■ 2(lu^,^ - 1] = 2mC/3 - fci , (24) 
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